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We discuss the generation and evolution of density perturbations during the large curvature phase 
of string cosmology. We find that perturbations in the scalar components of the metric evolve with 
cosmic time as exp("fH s t) where H„ is the Hubble constant during the string phase and 7 is a 
constant which is determined by an algebraic equation involving all orders in the string coupling a' . 
The seed for the perturbations can be provided by massive string modes. 



I. INTRODUCTION 



■ A very important issue in any cosmological model is the generation of density perturbations which evolved into the 
presently observed large scale structures of the Universe. At present, two main mechanisms have been studied. One 
is the amplification of zero point fluctuations at the transition between two cosmological phases, as in the transition 
from an inflationary phase to a standard radiation dominated phase; and the second possibility relies on topological 

■ defects produced at a GUT phase transition (see e.g. |l|,||). 

00 ' In the last few years, a cosmological model based on the low energy effective action of string theory has been 
proposed by Gasperini and Veneziano The model addresses the kinematical problems of standard cosmology 

thanks to a superinflationary 'pre-big-bang' phase (see however ref. Q for critical remarks) and to a subsequent Dc- 
Sitter like inflationary stage, and has a number of interesting phenomenological consequences ]^,|7j-4TT|] . The purpose 
of this paper is to examine the issue of the generation and evolution of density perturbations in this model. 

The analysis of refs. 12 indicates that the mechanism of amplification of vacuum fluctuations is not effective 



,^_< in string cosmology at very large wavelengths and it is probably unable to produce density fluctuations at the level 
5p/p ~ 1CU 5 as required for COBE anisotropy and for galaxies formation. Of course, it is still possible that density 
^Jy perturbations are generated by cosmic strings or other topological defects at the GUT phase transition, or in general 
• • . after the end of the high curvature phase of the model. In this paper, however, we explore the possibility that a 
. £h ' genuinely 'stringy' effect in the high curvature regime is responsible for the generation of density fluctuations. 

The paper is organized as follows. In sect. 2 we briefly recall the main features of the cosmological model. In sect. 3 
we compute the evolution of scalar perturbations during the high curvature regime in a simplified setting. (The details 
of the full analysis are relegated in the Appendix.) In sect. 4 we discuss a possible mechanism for the generation of 
the seed for scalar perturbation. Sect. 5 contains the conclusions. 



II. A SHORT SUMMARY OF STRING COSMOLOGY 



In this section we present the main features of the model and we establish the notations. For a more detailed 
presentation we refer the reader to refs. The model is obtained by considering the Neveu-Schwarz bosonic 

sector of string theory. The low-energy effective action is built from the massless states of the string: the graviton 
g^v, the dilaton cf> and the antisymmetric field and, in the string frame, it reads 

fi + (v*) 2 -^ff 2 l , (1) 
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where H^ vp — d^B vp ^ , R is the Ricci scalar]^ and A s is the string length. This action receives two kinds of correc- 
tions ||. The analogue of quantum field theory loop corrections, which are governed by e^, and genuinely stringy 
corrections which are controlled by the string constant a' (related to the string length A s and to the string tension \i 
by v^ 7 ^ A s = 1/^/Wp). 

In the pre-big-bang scenario the evolution starts from the string perturbative vacuum at low curvature and <fi — > — oo, 
where both kind of corrections are neglegible and physics is well described by the effective action ([!]), and it follows a 
phase of superinflationary evolution, which ends in a singularity at a finite value of cosmic time. However, when the 
curvature reaches a value of order one in string units, a' corrections to eq. ([!]), e.g. terms a' ' R 2 ^ vpa , become important. 
As shown in jlj|, at first order in a 1 they do indeed stop the growth of the curvature: the effective action at order a' 
receives a correction which depends on the Gauss-Bonnet term Rq B fl4[ and after a field redefinition at order a' it 
can be written in the form 



6 " 2A? ' " 



4, 



-ge 



R +(S/cf > f--(Rl B - (V^) 4 ) 



(2) 



where we have reabsorbed into a' a factor k = 1, 1/2 for bosonic and heterotic string respectively and we have assumed 
that the torsion background is trivial. 

The action (||) has a solution of the equations of motion of the form H = const. = H s , <j) = const. = c, where H is 
the Hubble constant, H Sl c are numerical constants of order 1/A S , and the dot is the derivative with respect to cosmic 
time. This solution acts as a late time attractor of the lowest order pre-big-bang solution obtained from eq. ([!]). Of 
course, when terms 0(a') are crucial, also terms 0(a' 2 ), etc. become crucial. Therefore we are in a truly stringy 
regime. The general structure of the equations of motions at all orders in a' is discussed in ref. [jl3||, where it is found 
that a solution of the form H = const. = H Sl <f> = const. = c, persists at all orders in a' if a set of two algebraic 
equations in the two unknown c,H Sl involving all orders in a', has real solutions. In the following, we assume that 
this is indeed the case. 

Finally, this 'string phase' should be matched to the standard radiation dominated phase. The problems connected 
with implementing this transition ('graceful exit') have been investigated recently by various authors pa]. 



III. EVOLUTION OF DENSITY PERTURBATIONS DURING THE STRING PHASE 

Metric perturbations around a FRW background can be classified according to their properties under spatial rota- 
tions. Scalar perturbations, which couple to density fluctuations, can be written in the longitudinal gauge in terms 
of two scalar functions ip% , ip2 as |l6| 

ds 2 = a 2 {j 1 ) {(1 + 2V>i)<V - (1 - 2V> 2 )dx 2 } , (3) 

where a is the scale factor, r\ is conformal time, and we are considering a spatially flat 4-dimensional FRW background. 
In this section, to illustrate our results, we limit ourselves to fluctuations with tpi — —ip2, i-e. to metrics which can 
be written in terms of a single function h as 

9ii.v{r},-x) = [l + /i(»7,x)] 3^(17) (4) 

where — a 2 ^)^^ is the background metric. The fact that these fluctuations are simply fluctuations in the scale 
factor allows a substantial simplifications of the relevant equations and allows us to present equations which are not 
too lengthy. As it will be clear below, we are only interested in general properties of the solutions, and we will show 
in the Appendix that these properties still hold if we consider the general case with ipi and ip2 independent. 

We also expand the dilaton as <f>(r),x) = <fto(r/) + <p(r),x). We make these substitutions into eq. (0) and expand 
to second order in h, (p. The condition that the linear terms vanish give the dynamical equations of motion for the 
background, 



R - (V0 O ) 2 + 2V 2 O - a' 



J (R 2 GB + 3(V0 O ) 4 ) - (V0 o ) 2 V 2 0o - 2V^oV^ V M V^ 



(5) 



*We use the sign conventions rj^u = (+, — , — , — ) and i?t pCT = d p T^ — . . . . 
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R - 2(V0 O ) 2 + 3V 2 0o - a'G^V^oV^o - V M V„<fo) = . (6) 

(The overbar denotes quantities constructed with g^, and G^ v is the Einstein tensor. Indices are raised and lowered 
with the background metric, and covariant derivatives V M are with respect to the background metric. The first 
equation is obtained with a variation with respect to <p and the second with respect to h.) Furthermore, the variation 
with respect to the lapse function gives a constraint on the initial values. 

For our homogeneous background, g^ v = a 2 ^)^,,, cf>o = 4>o(t]) the dynamical equations of motion, eqs. (pLm), become 



24> - % + 6^0- ~ &- 2 Q- a' (e^ + hi - tyffi - 3 2 <U = 0, (7) 
a a z a [as 2 4 a J 



30o - 20g + 90o- - S% - 6- - a' (-60 o ~~ + ~ 3< ^ ~ \ = 0, (M 

a a z a a a er a z a J 



while the constraint equation is 

, a 2 : o _ •, fl / f _ : a a 3 • 



6^-0 2 +60 o - + a' 60 o --t^ =0. (9) 
a z a L a 4 J 

Eqs. have the solution H — H s , <j) — c which, in terms of conformal time, reads 

ofo) = --!-, o (r?) = -^-log(-if s 7?) fa<0). (10) 

Notice that the ansatz (|o|) reduces eqs. (QJIJ)]) to 3 equations in 2 unknown c,H s ] remarkably, on this ansatz the 3 
equations are not independent and admit a common solution. This goes through at all orders in a' due to general 



properties of the beta functions of the underlying sigma model 13 1. 

Expanding the action at second order in </?, h and then varying it, or linearizing eqs. (^,^), we get the equations of 
motion for the fluctuations, which read 

2V^ V jLl <£> - 2VV + 3V 2 /i - 2VV V M /i -h(R- (V0 O ) 2 + 2V 2 O ) + 

-a' [-G^VpVvh + V m¥ > (4V^VT^o - 3V0o(V0o) 2 + 2V^ O V 2 O ) + (11) 

+v /1 v„v (2V0o V^o + g^i^M 2 )] = o , 

and 

-4V<foV M <p + 3V 2 p - 3V 2 /i + 3V"<£ V M /i + 

-a' [(V„V„/i) + VV0O - V^oV^o - ff""(V% - (V0 O ) 2 )) + (12) 
+(V^)(^i? - R^V V 4> ) + hG^iV^V^o - V^ VM) + 2V^G^V^ - G^V^V^] = . 

(In the above expressions we have not yet used the specific form of the background.) The linearization of the constraint 
equation is of course satisfied with the initial conditions h = tp = 0. 

Although formally correct, eqs. ( ^l| , ^2| ) suffer from a serious ambiguity. Consider for instance the coefficient of the 
term a'h in eq. (|l2]), i.e. '(V ' v 4>q — V M 0oV„0o)- We are free to add to it a term which vanishes because of the 
equations of motions of the background, e.g. an arbitrary constant times the left-hand side of eq. (||). The latter, 
however, is made of two separately non zero terms, one of order zero in a' and one of order a' . In a formal expansion 
in powers of a' of the coefficient of h in eq. Jl2]), part of the term that we have added would contribute to order a' and 
part to order a' 2 ; at the order at which we are working, this second part would be discarded. Thus, the coefficient of 
h in eq. (jlj), at order a', could be equally well written as 

- a' [G^(V M V„0o - V^oV^o) +a(R- (V0 O ) 2 + 2V 2 O )] , (13) 



with Ci an arbitrary constant. The same holds for all other operators in eqs. ( pTp^ ); for instance the coefficent of 
the term V M V„/i can be modified adding a term proportional to g^ v times the left-hand side of eq. (^). Therefore, a 
truncation of the coefficients at a finite order in a' gives ambiguous results. 

One might try to argue that there is at least a 'natural' prescription for the coefficients at any finite order which 
consists of taking eqs. (0,03) as they are, without adding terms proportional to the equations of motion. However, 



3 



even this is not true. Consider for instance the term (V 2 0o — (V0o) 2 ) which appears at order a' in eq. 
insert the explicit expression for fa, eq. (EOI), this is equal to 3cH s — c 2 . However, if we combine eqs. 




we find that (V 2 0o — i^fa) 2 ) is 0{a!) and, since this term appears in eq. (12) already multiplied by a', it gives a 
contribution 0(a' 2 ). Clearly, there is no natural choice between these two ways of treating this term. 

What we learn, therefore, is that a truncation at any finite order in a' gives completely arbitrary results. Only the 
full answer, including all orders in oi ', is meaningful. This is a consequence of the fact that the solution around which 
we are expanding emerges only after a' corrections are taken into account, and it is not a solution of the lowest order 
effective action, eq. (Q). So, the background solution mixes terms of different order in a' . 

We can however still use eqs. (El 12) in order to learn something on the structure of the equations at all orders. 



Using the explicit form of the background, eq. (10), we find that eqs. (|11|,|12|) can be recast in the general form 



rj 2 ( o «0 _ 40^) + afh + a ( S | + f _ 6 «v V ) + >?n% = , (14) 

where V is the fiat space spatial gradient and the index i = 1, 2 enumerates the two equations obtained from eqs. jTl| ) 
and ([12]), respectively. The coefficients a\ , etc. are independent of 77 and are formally given as a power expansion 
in a' . From eqs. ( |Tl|]l^ ) we formally get for instance = —3H 2 + a'(3H^ — cHg), etc., but as explained before, 
a truncation at finite order in a' is actually meaningless. Note that this would be true even if, for some numerical 
accident, H s , c were much smaller than one in units a' = 1. 

Let us introduce the modes hk(rj) performing the Fourier transform of h(r),x) with respect to x, and similarly for 



ip. Consider in eq. (14) the term aid 2 h/di] 2 — a2V 2 h, which can be rewritten as a2[{l/c 2 )d 2 hk/drj 2 + k 2 hk] where 
c.s = {0-2! a-\) x l 2 represents, physically the speed of the perturbation and, on general grounds, c s < 1. For wavelengths 
with r]kc s < 1 the gradient term k 2 hk is small compared to (1/ 'c 2 )d 2 hk/ 'drj 2 ~ hk/{c 2 rj 2 ). Of course, this is just the 
condition that the physical wavelength of the perturbation is greater than the Jeans length. On the other hand, if 
c s < 1, there is a window of values of the wavelength for which at the same time the perturbation is still within the 
horizon, krj > 1. The latter condition is important because for super- horizon-sized modes it becomes relevant the fact 
that the scalar perturbations ipi, ip2 are not exactly gauge- invariant. Indeed, ipi, 1^2 are the Bardeen variables written 
in the longitudinal gauge [ pT|]l(| , and are gauge-invariant only at first order for small coordinate transformations. 
The effect of non gauge invariance becomes important for physical wavelength larger than the horizon, A p h ys > H , 
where A p h ys = 0(77) A. In a DeSitter background this means krj < 1 where k = 1/A is the comoving wavenumber. For 
krj < 1 a growing perturbation can be a gauge artefact with no physical meaning. 

If instead 1 < krj < l/c s there is a range of wavelengths for which all gradient terms can be neglected but the 
perturbations are still within the horizon; then we see that equation (^) for hk,ifik depends on 77 and d/drj only 
through the combination rjd/drj. This was not obvious a priori. On dimensional grounds, 77 could also have appeared 

(i) 

in combinations as r]H s ,rjc, and therefore the coefficients a\ , etc., could have been functions of rj, while instead they 
only depend on the constants H s , c, a 1 . This result can be traced to the specific form of the De Sitter background and 
of the dilaton, eq. (|ic|). To understand this point, observe for instance that in this background, using conformal time, 
the Ricci tensor is given by R^ v = —ZH^r^rj^, while g^ v = H 2 rj 2 rj^ v , d^fa — g^ LV d v fa = —cH s rj5Q and similarly 
V^V^^o ~ V 2 - So, in a term like V M V„/i in eqs. ( |ll| , p^| ), independently of whether we use W" v , g^ y or derivatives with 
respect to fa to saturate the indices ji, v, we get a factor of 77 for each index to be saturated. One can easily check 
that the argument goes through for all operators one can construct, and therefore it holds at all orders in a' . For 
instance, at higher orders we will have terms with more than two derivatives, e.g. V^V^ \/ p V a h. Whatever operator 
we use to saturate the four indices, we get a term ~ rj A d A hj drj 4 . 

The existence of a range of values of k for which the perturbation is within the horizon, krj > 1, but still the gradient 
term can be neglected, is a general properties due basically to the fact that the speed of the scalar perturbation ('speed 
of sound') is smaller than the speed of light, and we expect that it goes through at all orders in a' , when higher order 
spatial and temporal derivatives come into play. If this is the case then, for such wavelengths, the full equations 
governing the perturbations, at all orders in a' , have the general form 

Fi(v^ ] )hk(v)+G i (v^)Mv)=0, (15) 

where i — 1,2 enumerate the equation and Fi, Gi are unknown functions. The non-trivial content of eq. (03) is that 
they depend on 77 only through the combination rjd/drj. Looking for solutions of the form hk ~ 77 7 , ip^ ~ 77^ reduces 
eqs. ([15]) to a set of two algebraic equations for the two unknown /3, 7: 

Fi( 7 )+Gi(/3) = 0, F 2 (-f) + G 2 ((3) = 0. (16) 
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The situation is particularly interesting if the above equations have a real and negative solution for 7. In fact, since 
in DeSitter space the relation between conformal and cosmic time is —H s rj = exp(— H s t), we get an exponentially 
growing mode 

h k ~ exp{\i\H a t} . (17) 

The maximum amplification which can be obtained is related to the duration of the string phase. Denoting with rj s the 
value of conformal time at which the string phase begins and with 771 the value at which it ends, (rj s , rji < 0,\r] s \ > \r]i\) 
the maximum amplification factor is (77^ / 771 ) . The parameter (rj s /r)i) (equal to fi/f s in the notation of ref. ||) is a 
free parameter of the model, and can be very large. In fact, the most interesting phenomenological situation from the 
point of view of the production of relic gravitational waves is realized when this number is at least of order 10 8 

We see therefore that string cosmology has a possible built-in mechanism which allows to amplify in a substantial 
way very tiny initial inhomogeneities. 

Of course, the above result only holds in the linear regime. However, this is sufficient for our purpose, which is to 
find a mechanism which generates density fluctuations at the level (Sp/p)k ~ hk ~ 10 -5 at wavelength of the order 
of the present Hubble radius or at the scale of galaxies formation. 



IV. THE SEED FOR PERTURBATIONS 



In the previous section we found that, depending on the form of some unknown algebraic functions Fi,Gi, which 
involve a full (all-orders in a') computation, we can have a very large enhancement of initial density fluctuations. The 
next step is to identify a possible candidate mechanism which generates a seed which is then amplified. 

In the dilaton dominated regime where the action (Q) gives a good description of physics, there is no mechanism 
which can generate inhomogeneities. In fact, as recently shown by Veneziano even starting with generic inhomo- 
genehous initial conditions the model evolves into a highly homogeneous spacetime. Therefore, we have to look again 
into the string phase in order to find a seed for density perturbations. 

Until now we have limited our attention to the massless modes of the string. However, massive string modes are 
associated with a' corrections to the action fH^ , po|| , and must be included for consistency once we include corrections 
of the type a' R 2 ^ vpa . In fact the complete effective action of the string has the general form [pof 

00 N n 

S= f dVJ^aT E Og>(dX)Bg\x), (18) 

J n =0 i„=l 

where the string coordinates have been rescaled, A M — > \fa'X tl , Bi n are the background field at a given massive 
level (rescaled so that they are dimensionless) and the corresponding operators; N n is the number of irreducible 
operators for given n. The power of a' is uniquely connected with the mass level. At least working order by order in 
the mass level [^(J , the condition of quantum Weyl invariance of the action ( |l8| ) induces a' corrections to the equations 
of motion derived from eq. (Q). Since we have seen that a' corrections have a crucial effect in the string phase, it 
is clear that in this phase it is not appropriate to restrict the attention to the graviton-dilaton-antisimmetric tensor 
field sector. 

Of course, a description of the full stringy regime is very difficult to obtain. However, a number of general properties 
are well understood from the study of strings at very large temperature and densities E| . The main results are as 
follows. At low temperatures and densities the canonical ensemble provides a good description of a gas of strings. 
As T approaches the Hagedorn temperature Tjj, the canonical ensemble breaks down and one has to resort to the 
more fundamental microcanonical ensemble. As T — > T^j the canonical energy density tends to a finite value p c (in 
dimensions d > 4). The most intriguing result is that if we increase the energy density beyond this critical value, all 
the excess energy goes into a single highly excited string, rather than raising the temperature of the string ensemble. 

In the high curvature phase of string cosmology the parameter which controls the value of the curvature is H s . If 
H s exceeds a critical value H c , the energy density of the gravitational field exceeds the critical energy density p c and 
the excess energy goes into a very excited string. More precisely, we should expect the formation of one very excited 
string in each horizon volume, since regions separated by more than a horizon distance behave independently.^ This 



Here we are using the common abuse of language of calling H 1 the horizon distance. Actually, in string cosmology the 
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means that, if H s > H c , the favored thermodynamical configuration is highly inhomogeneous, and the massive string 
modes can act as the seed for density perturbations. 

To understand more in detail the effect of the highly excited strings, let us denote by 0^ their energy momentum 
tensor. It acts effectively as an external source, and the equation of motion of the perturbations obtained with a 
variation with respect to g^ or, in our case, with respect to h (which corresponds to i = 2 in eq. Jl5|)) becomes, at 
all orders in a', 

F 2 ( V ^-)h k ( v ) + G 3 (r,£-)<p k (ji) = A S V9 , (19) 
arj di] 

where = g^Q^, while the equation with i = 1 remains homogeneous. 

As already mentioned, eqs. (||J6|), and therefore their linearization, eqs. ( |ll| , [T2| ), do not exaust the content of the 
classical theory. Another equation is obtained varying the action with respect to the lapse function N, which is then 
set equal to one with a choice of gauge. This equation, of course, is a constraint on the initial data, which is conserved 
by the dynamical equations of motions. The linearization of the constraint equations gives as the initial conditions for 
the perturbations h = p — (and similarly all the time derivatives of h, ip can be set to zero up to one order less than 
the higher order derivative appearing in the dynamical equations). Since the energy momentum tensor is coupled to 
the lapse function, the constraint equation in the presence of massive string modes has the general form C(h, ip) = Q 
where C is a function of h, p> and of their derivatives which vanishes when all its arguments are zero. 

Therefore, we see that the 'external source' has a two-fold effect. First, it forces h to have a non-zero initial 
value, through the constraint equation. And second, the evolution of the perturbation will be the sum of a particular 
solution of the inhomogeneous equation Jl9| ) and of the general solution of the associated homogenehous equation, 
which includes (if 7 < 0) the exponentially growing mode ([l^). These two effects, which are quite general, correspond 



to what have been termed as 'initial compensation' and 'subsequent compensation' in ref. \22\. 

It is quite difficult to estimate the form of the spectrum of density perturbations, since it depends on the details of 
the string phase, and we leave the problem open for further work. Even more difficult is to estimate the numerical 
value of (Sp/ p)k, since it involves in an essential way the numerical values of rj s , 771,7 on which we have no clue. 

V. CONCLUSIONS 

The problem of the generation and evolution of density perturbations in string cosmology appears to be a very 
difficult one, since it fully involves the details of the high curvature stringy regime, for which we do not have an 
adequate theoretical understanding. 

Still, a number of general considerations appear to be possible. It turns out that in a number of situations some 
dynamical questions concerning the string phase and which therefore in principle are very complicated, can actually 
be reduced to the issue of the existence of a solution to a certain algebraic equation involving all orders in a' . 

An example of this phenomenon has been found in ref. [T^ ] , where the issue of the existence of a solution of the 
equations of motion which forbids the lowest order pre-big bang solutions from running into a singularity has been 
shown to reduce, to all orders in a' , to the existence of a real solution of an algebraic equation. In this paper we 
have found a similar characterization for the existence of exponentially growing scalar perturbations during the string 
phase. The basic mechanism is that, for wavelengths larger than the Jeans length, in a DeSitter background with a 
dilaton growing linearly with cosmic time, the equations of motion of the perturbation have solutions of the form 7y 7 
and, in DeSitter space, 77 ~ exp(— H s t); this results, for 7 < 0, in an exponentially growing mode. 

From the physical point of view, both the exponential growth of the perturbation, and the generation of seeds from 
highly excited strings appear under rather general circumstances and it is therefore possible that string cosmology 
has its own, genuinely stringy mechanism for generating large scale density perturbations. 

As already stressed in |Q , an all order solution should correspond to an exact conformal field theory (see p3j for 
work along these lines) and our algebraic conditions could have a natural counterpart in the conformal field theory 
approach. 

We thank Maurizio Gasperini and Gabriele Veneziano for useful comments. 



particle horizon is infinite H], but still H 1 gives the scale of the distance between points which formally recede from each 
other at the speed of light. 
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APPENDIX A: 



In this Appendix we examine the equations governing the evolution of the perturbations in the general case with 
V>i and ip2 independent. Technically, we lose the advantage of having a metric conformal to Minkowski metric, and as 
a consequence most of our equations become very lengthy. However we are only interested in the general properties 
of these equations, and in particular in the fact that they only depend on rjd/dr] once spatial gradients have been 
neglected. 

The equations of motion for the background are of course unchanged. We now have three dynamical equations 
of motion for the fluctuations, which can be obtained by expanding the action at second order and performing the 
variations with respect to (p t ipi,ip2- The variation w.r.t. ip gives 



2V /1 <pV 



2VV + 2V M (V>i + V2)VVo + 4V (V>i + $0V% - 4V M ViV 1 </> - 6V 2 V>2 + 2V fc V fc (V>i + V2) + 



+(V>i + ty 2 )R% - (^1 - + (V>2 - V-i)(2V 2 o - V M V^ ) + 2(V>i + V 2 )(2V V> - Vo^V%) + 



— {SV^G^ - R 2 g [ b ] + 3R 2 GB {^ + V2) + 4(V0 O ) 2 (V 2 ^ - 3V^ V» + 16V„VVoVVoV^ - 

-8V 2 </) V^V^ + 8V M VVV^oV^o + (V>i + V> 2 ) [-36V V%V 6>V% + 9(V0 O ) 4 + 
-4ViVVoV ^oVVo] - 12V (^i +^2)(VVoV o 0oVVo)} = 0, 



(Al) 



where repetead latin indices like k. . . mean summation over the three spatial coordinates and the index denotes 
the time coordinate. The variation w.r.t. Vi gives 

2V l VV - 2V^V^ + 3(V>i + ^2)R° - (V>i + ^2)Rl + 

+(Vi + V2) [(V0 O ) 2 + 2V o 0oV% - 2V fc vVo] + 4V fc V fc V 2 - 2V^ V^V2 - 4V ^oV°V>2 + 
-j{bp5& B + 9R 2 GB {^2 - V>i) + 40G^V M V^2 + 32V 2 $ 2 R - 8V0V2VV0-R + 8(V>i + V^X^Oi + H) + 
+4(V>i + V2) (i?f;) 2 - 8(V>x + V2) (R° k )) 2 + 8(Vi + V2) (R k *) 2 - 8 (ViV^Vi + V2)) - 8(V>x + V2) (i?S) 2 + 
-16(Vi + V2) (-R}) 2 - 32(^i + 1P2) (R°) 2 + 8 (ViV J '(V»i + ^2)) i?} + 4(Vi + ^2)^ + 



+4 (ViV*^! + V2)) (£g ~ Rt) + V>1 8R^(R^i + Ri) + 4(^) 2 - 8(i?«) 2 + 8(^j) 2 



-8 JZ 



fl? V,-V 



2 



+ 



(A2) 



8(i?°) 2 - 16(R{) 2 - 32(i?j,) 2 + 4RRj - 4 (V fc V fc 0o) (i?o - R 
+8(5^i - V2) (-R°i) Z + flg? [16VjVVi ~ 8V,VV 2 + 16V J V>o(2V'2 - V'l) - 32^^-] + 

+32(4^ + V2) (R%) 2 + R° [-24V 2 V 2 + 40V fe V fe V2 - 16V fe V fc Vi + 36V^2V> - 8^(5^1 - V> 2 ) + 
+16^iV fc V fc o - 8^|(4^i + V2) - 24V> 2 V fe vVo] + R) [-56V;V J 'V 2 - 16V 4 V>i + 16(Vi - 2V> 2 )V i V i V>o] + 

[28V 2 V- 2 + I6V0V2VV0 - 24V-1-R* + BthVjV'fo] + R [-48V 2 V> 2 + BViVtyi - 8(V>i - V^ViVVo + 24Vi#] + 
-8V 4 V J V,V l Vi + 8 (ViV J '(^i - V2)) (V,V>o) - 16 (V,-V J '^) (V;VVo) + 24 (V 2 V> 2 ) (V;V>o) + 
-24V 2 V l V V°V 2 + 24i?^V A1 V> 2 - 24 (V"Va) iffivfo + 

-3(V o 0oVVo) 2 (5V'i + 3^2) + 12(V o ^V%)(V o 0oV%) - 3^(V o 0oV%) 2 } = 0. 

To simplify the equations, in the terms between curly brackets we have used the fact that cf>o does not depend on 
spatial coordinates, and that during the stringy phase every component of Riemann and Ricci tensor with an equal 
number of upper and lower indices is costant in time as well as in space. 

The equation that one would obtain by varying the second order action w.r.t. V2 can be obtained more simply 
noting that subtracting it from eq. (A2) one gets equation (|l2]). 

Having all the equations, we can now directly check by inspection that the property of homogeneity is indeed 
preserved. If we neglect all spatial gradients, the equations (Al), (A2) and jl2] ) have the general form 

Mv^M + B^^-^kiv) + Ci(v^)Mv) = 0, (A3) 



where now i — 1,2,3, and therefore we can look for solutions V'l ~ ?? 7l ,V'2 ~ V 12 ^ ~ V 13 \ which reduce eqs. (A3) to a 
set of 3 algebraic equations for 3 unknowns 71,72,/?. The general situation is therefore completely analogous to that 
discussed in the text. 
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